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Abstract 

In this paper we study to what extent the canonical equivalence and 
the identity of the geometric phases of dissipative and conservative linear 
oscillators, established in a preceeding paper, can be generalized to nonlinear 
ones. Considering first the 1-D quartic generalized oscillator we determine, by 
means of a perturbative time dependent technic of reduction to normal forms, 
the canonical transformations which lead to the adiabatic invariant of the 
system and to the first order non linear correction to its Hannay angle. Then, 
applying the same transformations to the 1-D quartic damped oscillator we 
show that this oscillator is canonically equivalent to the linear generalized 
harmonic oscillator for finite values of the damping parameter (which implies 
no correction to the linear Hannay angle) whereas, in an appropriate weak 
damping limit, it becomes equivalent to the quartic generalized oscillator 
(which implies a non linear correction to this angle) . 
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INTRODUCTION The Hannay's angle [1] (classical counterpart of the Berry's 

geometric phase [2]), originally associated with the adiabatic evolution of 
classical hamiltonian systems, has been recently extended to a large class of 
dynamical equations corresponding to dissipative systems: non linear equa- 
tions with limit cycles [3] or with more general internal symmetries [4] , equa- 
tions describing the dynamics of the laser [5], etc... In this context we have 
shown in a recent paper [6] that the simplest dissipative system, namely the 
damped harmonic oscillator specified by the time dependent Lagrangian 

C{qA^) = \e 2 S tx{s)ds {q 2 -ulq 2 ) (1) 

is canonically equivalent, even for time dependent parameters (X,cu ) = fl, 
to the generalized harmonic oscillator, a conservative system specified by the 
Hamiltonian 

H(P,Q,[l) = ^- + \PQ + ^-Q 2 . (2) 

(The generating function of the canonical transformation is F(q, P, t) = 

qPel x(s)ds ). As a consequence the Hannay's angles of the two systems 
are identical. Their expression can be simply recovered, like in [7], using the 
transformed variable 

Q= dF{ %p' t) =qefM* (3) 
which brings the Lagrangian (1) to the form 

HQ, Q, = \(Q 2 - 2X QQ ~ " 2 Q 2 ) (^ 2 = ^ - A2 ) ( 4 ) 

which also reads 

L(Q, Q, ft) = \{Q 2 - (oo 2 - X)Q 2 ) - \j^Q 2 ) ■ (5) 

In this later expression the quantity (cu 2 — A) is the square of the instantaneous 
frequency of the system. In the adiabatic limit where the parameters \x are 
slowly time-varying functions fl(et) (with — C 1) one can make a first order 



3 



expansion of this instantaneous frequency with respect to the small adiabatic 
parameter e. One gets the well known result [1,2] 

a=«-A. (6 ) 

where the 'dynamical' part u of the time derivative of the phase of the 

oscillator appears to be corrected by an adiabatic contribution — — the 

integral of which in the parameters space is the 'geometrical' Hannay's angle. 

The main purpose of this paper is to study to what extent the above re- 
sults, concerning the canonical equivalence and the identity of the geometric 
phases of dissipative and conservative linear oscillators, can be generalized to 
nonlinear ones at least in the first order approximation of the perturbation 
theory. In the following we restrict ourselves to the quartic generalized oscil- 
lator and compare it with the quartic damped oscillator. We do not consider 
cubic terms because, at this order, they are non resonant i.e. without effect 
on the phase equation [8]. 

The first section is devoted to the quartic generalized oscillator the Hamil- 
tonian of which is deduced from (2) by the addition of a term proportional 
to Q 4 . In order to calculate the geometric phase of this system we extend 
the method of reduction to normal forms to the case where the parameters 
vary slowy with time. We are then able to solve perturbatively the Hamilton 
equations and to find the appropriate canonical transformations under the 
two assumptions of weak nonlinearity and of adiabaticity of the variation 
of the parameters. In particular we determine the adiabatic invariant / of 
the system and we find the expression of the first order nonlinear correc- 
tion to the geometric Hannay part of the angle 0. The reason for choosing 
the normal forms technic rather than the averaging method is that, besides 
the fact that it can in principle be developed to any order of perturbation 
(in the adiabatic regime), it is a standard approach for the study of nonlin- 
ear dissipative systems; it will allow us to deduce the appropriate canonical 
transformations in the dissipative case from the above ones. 

The second section is devoted to the study of the damped quartic oscilla- 
tor the Lagrangian of which is deduced from (1) by the addition of a term pro- 
portional to g 4 . Using the time-dependent canonical transformations found 
in section 1 and keeping the same hypotheses of weak nonlinearity and of 
adiabaticity we obtain the following results. For finite values of the damping 
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parameter A, i.e. when resonance does not exist, the quartic term is without 
effect on the geometric part of the angle. In that case the quartic damped 
oscillator can be shown to be canonically equivalent to the linear generalized 
harmonic oscillator. However, there also exists a weak damping limit char- 
acterized by a magnitude of A going to zero with the adiabatic parameter 
in such a way that the resonance phenomenon resurges in the limit. In this 
limit the quartic term contributes and one recovers the Hannay's angle of the 
quartic generalized oscillator determined in the previous section. This later 
result generalizes to nonlinear systems the result established in [6] for linear 
oscillators. 



1. GENERALIZED QUARTIC OSCILLATOR The simplest nonlinear ex- 
tension of the generalized harmonic oscillator leading to a resonant term in 
the equation of motion is obtained by adding a quartic term to the Hamil- 
tonian (2) which thus becomes: 

H G (P, Q, # = ^ + XPQ + ^Q 2 + V -Q\ (7) 

The Hamilton equations for Q and P read 

Q = XQ + P , P = -XP - ulQ - uQ 3 . (8) 

In order to solve these nonlinear coupled equations it is convenient to intro- 
duce, in place of P and Q, the complex variable z, and its complex conjugate 
z, defined by 

:7r Q--(AQ + P)l, (c 2 = ^ 2 -A 2 ). (9) 

Equivalently Q and P read: 

1, . ^ ius — X iu> + X 
Q = -f=(z + z), P = ^=z i=-z . 10 



For time dependent parameters (A, ujq, v) = jl, the Hamilton equations for Q 
and P lead to the following nonlinear equation for z: 

iv , iX Co — iX 
z = iujzA -(z + zY z. (11) 
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(Note that the three parameters A, u and v do not play the same role: the 
time derivative z> of the parameter associated with the non linear quartic term 
does not appear in (11) in contradistinction with A and to.) This equation 
can be solved perturbatively using its canonical reduction to normal form 

[9] if one assumes that the system is weakly nonlinear (— Q 2 <C 1) and that 

° 

the parameters vary adiabatically (/2 is a slowly time varying function ft(et) 
with — < 1). Under these conditions, let us introduce the near to identity 
transformation 

z = u + 5u (12) 
in order to eliminate the nonresonant term proportional to u into the equation 

for the new variable u. For 5 = - ' ' (8 small, of order e) the coefficient of 

Auj 2 

u cancels and the equation for u reads: 

u = i(w- ^)u+^(u + u) 2 [(l + 8 + 35)u+(l + 4:5)u\ . (13) 

The first term in (13) already exhibits the Hannay's angle of the linear gen- 
eralized oscillator. In order to obtain the non linear correction one must 
introduce a second (near to identity) change of variable 

u = v + av 3 + (3vv 2 + >yv 3 . (14) 

The requirement that the equation for v no longer contains the nonreso- 
nant cubic terms present in (13) leads to differential equations for the time 
dependent coefficients a, (3 and 7 of the transformation. These differential 
equations are explicitly written and their solutions discussed in the appendix. 
The resulting equation for v, valid up to the first order in e and in the weak 
nonlinearity parameter's u, is: 

i = ^-l )v + ^ {1 + ^ )v2Ti (15) 

Then, setting v = Ae te , the equations for the amplitude A and for the angle 
6 read: 

A = 0, (16) 

<d=u(l + ^-A 2 )-±(l-^-A 2 ). (17) 
V 4w 3 > 2u V 2a; 3 > K J 
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Equation (16) shows that A is an adiabatic invariant related, as seen below, 
to the action I of the system 

I = A 2 . (18) 

The first term in the r.h.s. of equation (17) accounts for the well known 
result that the quartic term —Q 4 in (7), responsible of the presence of the 

resonant v 2 v terms in (15), induces a renormalization of the linear frequency 
uj which becomes 

n = u(l + p-A 2 ) . (19) 
4a; 3 

(The facts that the frequency does not explicitly appear in the expression 

of / and that the correction term — - A 2 to the linear frequency is different 

Auj 2 

from the expression found in reference [8] are due to the use of the amplitude 
A of the transformed variable v in (18) and (19) in place of the amplitude a 
of the original variable Q; in particular, the usual expression for Q in terms 

UJ 



of a, uj and v can be recovered from (19) noting that A — J — a.) The second 

V 2 

term — (—1 + -^-^A 2 ^j in the r.h.s. of equation (17), which exists only for 

time dependent parameters, is the geometric, non integrable, Hannay's part 
of O. Like the dynamical part Q it also contains a contribution from the 
nonlinear resonant term. However, as this geometrical part is defined up to 
a total time derivative one can write it, for v constant, under the same form 

as in the linear case ( — 7^~~) m terms of the renormalized frequency £1 and of 
a renormalized damping parameter A defined by 

A = A(l-^A 2 ). (20) 

T , , 1 ,. A A . 3v 9 . d /9u\ 9 \ . 

Indeed, the equality — — = — 1 A ) — - A is exact up to the 

4 3 2Q 2uj k 2uj 3 1 dt^8uj 4 > l 
first order in the weak nonlinear term proportional to vA 2 . Consequently, 

the equation (17) for the angle variable 9 takes a functional form, 
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identical to the one, O = to — — , obtained in the linear case. One can 

2uj' 

say that, for v constant, the effect of the weak quartic nonlinearity on the 
phase (the 'geometrical' as well as the'dynamical' parts) amounts to a simple 
renormalization of both uo and A. 

We now show that the dynamical variables (J, 0) and (P, Q) are related 
by a time dependent canonical transformation and that the equations (16- 
18) are the Hamilton equations for action-angle variables associated with the 
time dependent Hamiltonian 7^g(J, 0, t): 

H (/,e,t) = /o,(i-A) + / ^( 1 + ^) . (22) 

This can be verified by inspection of the set of successive transformations 
(P,Q) — > (—iz,z) — > (iu,u) — > (—iv,v) — > (I, 0), which will be also of 
importance in section 2. As concerns the first transformation (P, Q) — > 
(— £z, z) its generating function F(Q, z) is obtained by integration of the 

OF dF 
equations — — = P(Q, z) and — — = iz(Q, z) deduced from the differential 

identity (characterizing a canonical transformation) 

PdQ - H G dt = -izdz - H z G dt + dF . (23) 
Taking into account the relations (9-10) one gets 

F(Q, z) = -hiu + X)Q 2 + iV2uQz - \%z 2 . (24) 

The Hamiltonian H G for the new conjugate variables (—iz, z) is then obtained 

dF 

from the relation 7i G = H G + — — . Its expression reads 



nh = ^ + jb iz + ' )4 + ^ -^-h z+r>2 (25) 



and one can verify that the Hamilton equations 

&H Z G ^ &H G 



a = — (26) 



d(—iz) dz 

indeed coincide with the equation (11) for z and the corresponding one for 
z. 
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The transformation (—iz, z) — > (iu, u) associated with the generating 
function 

G(z, u) = -izu + -i5u 2 - -5z 2 (27) 

is also canonical. The expression of the Hamiltonian H G for the new conju- 

dG 

gate variables (iu, u) is obtained from the relation H G = Ti G + and, in 

this hamiltonian formulation of the reduction to normal form, 6 and 5 are de- 
termined by the requirement that Ti G no longer contains terms proportional 
to u 2 and u 2 . This is equivalent to the above requirement that the equation 
for u (resp. u) no longer contains non resonant term proportional to u (resp 

u) and it leads to the same value - for S. Then, Hr reads 

Alo 2 

n u G = (u-^)uu+-^(u + u) 4 + ^(u + u) 3 (5u + 5u) . (28) 
dG 

(In fact -7^- does not contribute to this expression of H G valid up to the first 

order in e since it only contains terms proportional to 5 and 5 and is thus of 
order e 2 .) 

In the same way, the transformation (iu,u) — > (—iv,v) associated with 
the generating function 

1 11 

K(u, v) = ivu + av 3 u + -(3vu z + ^T^ 4 — ^T^ 4 (29) 

is canonical. The values of a, (3 and 7, now determined by imposing that the 

dK 

Hamiltonian TL V G = TL G + -7^- for the new conjugate variables (— iv, v) does 

not contain quartic non-resonant terms, are those found in the appendix; the 
expression of H G is 

«& = <«-^+j£(l + >V. (30) 

Finally the transformation (—iv, v) — > (I, 0) associated with the gener- 
ating function 

M(v, 6) = --v 2 exp(-2iQ) (31) 
2 



is canonical and the Hamiltonian (30) transforms into (22) as announced. 
This expression (22) shows that / = A 2 is the action of the system and 
allows to determine the adiabatic invariant as a function of the energy. 2. 

DAMPED QUARTIC OSCILLATOR 

The damped quartic oscillator which we consider in this section is de- 
scribed by the equation: 

q + uj 2 q + 2Xq + uq 3 = 0. (32) 

It can be obtained from the generalized time- dependent Caldirola-Kanai La- 
grangian [10] 

C D (q,q,P) = \e 2 ^ s)ds {q 2 -^q 2 - \^). (33) 
In terms of the new variable 

Q = qef x ^ ds (34) 

the Lagrangian (33) reads 

Ld(Q, Q, fl) = \(Q 2 ~ - uj 2 Q 2 - ±i/QV 2 (35) 

and the corresponding time- dependent Hamiltonian H^(P, Q, p), where P = 

— ?- = Q — AO, takes a form : 
dQ 

h d (p, Q,fi = Y + Xp Q + y Q 2 + \Q^^ tx ^ ds (36) 

identical to the expression (7) for the Hamiltonian of the quartic generalized 
oscillator, except the crucial exponential factor in the quartic term. Then, 
making the same successive changes of variables as in the previous section, 
namely {P,Q) — > (—iz,z) — > (iu,u) — > (— iv, v), we obtain without any new 
calculation the following Hamiltonian for the conjugate variables (— iv, v): 

W D = {u- ±)vv+ p-Jl + - 2 )v 2 v 2 (37) 
Zlo our ur 
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At this point a remark is of order. In contradistinction with the gener- 
alized quartic oscillator the damped quartic oscillator does not exhibit the 
resonance phenomenon. The theory of normal forms teaches that it is possi- 
ble in that case to eliminate non linear terms. As a consequence one expects 
(as announced in [6]) that the Hannay's angle does no get any non linear 
contribution. Let us see how this comes in the formalism of canonical trans- 
formations. To this end we introduce the (near to identity) change of variable 

v — w + iaw 2 W. (38) 

The transformation (—iv, v) — > (iw,w) is canonical for a real. It corresponds 
to the generating function 

N(v,w) = —ivw crv 2 w 2 (39) 

dN 

and the transformed Hamiltonian Hp = H V D + -j— reads 

C6 

HI = (a; - ±)v*D - l -(a-^- 2 {l + A) e -2rAM*) w V • (40) 

The second term in the r.h.s. of (40) can be set equal to zero. Indeed, 
because of the presence of the exponential factor the corresponding value of 
a remains small (proportional to v) which guarantees that the transformation 
(38) is close to identity. Thus for A finite the quartic damped oscillator is 
canonically equivalent to the generalized harmonic oscillator and the phases 
of the two systems are identical. Note that such an elimination can not be 
done on the Hamiltonian (30) because the integration of a would lead to 
unbounded terms proportional to the time t. 

Let us now consider the case where the magnitude of the damping param- 
eter A is close to zero in such a way that the resonance phenomenon cannot 
be avoided. More precisely let A be of the form 

\(t) = e a \(e b t) (41) 

where a and b are positive numbers such that a > b > and a + b — 1. b 
positive ensures the validity of the adiabatic hypothesis for A and a + b = 1 
that A remains of order e, like Co. For such a behaviour of A the integral 
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/*A(s) ds is of order e a ~ b with (a — b) > 0. Then, replacing the exponential 
factor in (37) by the relevant terms of its expansion one gets the following 
expression (valid up to the first order in e) for the Hamiltonian TC V D : 

n v D = (u - ±)vv + |^(l-2 f\(s) ds + \)v 2 v 2 . (42) 
Zui OUT J co z 

For the same reasons as for the Hamiltonian (30) the non linear terms in (42) 
cannot be eliminated. It is interesting to note that the term proportional to 

A 3z/A 

the time derivative of the parameters, vv-\ jv 2 v 2 , is identical in (30) 

and (42). Therefore in this weak damping limit the Hannay's angle of the 

quartic damped oscillator is the same as the one of the quartic generalized 

harmonic oscillator. Due to the presence of the supplementary dynamical 
3v ft 

term — - — - / A(s) ds v 2 v 2 in (42) the two systems seem to be not canonically 

equivalent. However one can express the Hamiltonian (42) in terms of the 
renormalized parameter v = v{l — 2 J*A(s) ds) (a change which, in the degree 
of accuracy of our calculations, does not affect the geometrical nonlinear part 
of the angle) to get an expression 



Til = (u - -)vv+ — (1 + -)v 2 v 2 (43) 



A , 3i> . A 
2^ + 8^ (1 + ^ 
similar to (30). The damped quartic oscillator with parameters Uq, A and v 
is thus, in this weak damping limit and in the first order approximation of 
perturbation theory, canonically equivalent to the generalized quartic oscil- 
lator with parameters uj , A and v and both systems have identical Hannay's 
angles. These results are the generalization for nonlinear oscillators of the 
result established in [6] for linear ones. 

APPENDIX 

The nonlinear, near to identity, change of variable (14) transforms the 
equation (13) for u into the following equation for the new variable v: 

v = i(u-—)v+-^(l + 2(5 + 5))v 2 v-D 1 a v 3 - D 2 /3 vv 2 - D 3l v 3 (Al) 
v 2uj ' Auj 2 

where the differential operators D 1: D 2 and D 3 are such that 

D 1 a = a + 2i(u- — )a- ^(1 + 8 + 35) (A.2) 
v 2u J Au 2 
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D*P = P- 2<(w - — )0 -— 2 (l + 35 + 5) (A3) 

D 3 7 = 7 - 4*(cu - A) 7 _ £_(i + 45) (A4) 

For a solution of Dia = 0, /3 solution of = and 7 solution of D37 = 
the nonresonant cubic terms can be eliminated in (A.l). The three differential 
equations having the same structure, it is sufficient to study the behaviour of 

the solution of one of them, for example a. Since S = — — — , the differential 
equation for a reads 

Q + 2l ( w __) Q __( 1 + __) = 0. (A5) 

This equation which contains terms of order zero and of order one with 
respect to the small adiabatic parameter e can be solved by a perturbative 
method and its solution can be written under the form of an expansion with 
respect to e. Neglecting the terms of order e in (A. 5) (i.e. all the time 
derivatives) one finds the zero order approximate solution a = — Then 

putting a = «o + ol\ into (A. 5) one gets an equation for a\ which contains 
terms of order one and of order two with respect to e. Neglecting the terms 
proportional to e 2 in this equation, one obtains an expression (of order e) 
for a.\ and so on. One verifies that a, and thus also (3 and 7, are small 
for weak nonlinearity and in the adiabatic hypothesis (the leading terms are 
proportional to v and e). This validades the above calculations which suppose 
that the transformation is near to the identity. Finally the nonresonant cubic 
terms can indeed be eliminated in (A.l) which reduces to (15). 
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